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SUMMARY 

Current literature on the three-dimensional flow through compressor cas- 
cades deals with a row of rotor blades In Isolation. Since the distance 
between the rotor and stator Is usually 10 to 20 percent of the blade chord, 
the aerodynamic Interference between them has to be considered for a proper 
evaluation of the aerothermodynamlc performance of the stage. A unified 
approach to 'the aerodynamics of the Incompressible flow through a stage Is 
presented that uses the lifting-surface theory for a compressor cascade of 
arbitrary camber and thickness distribution. The effects of rotor-stator 
Interference are represented as a linear function of the rotor and stator flows 
separately. The loading distribution on the rotor and stator blades and the 
Interference factor are determined concurrently through a matrix Iteration 
process. 


INTRODUCTION 

The multistage axial compressor functions In such a way that each stage 
performs essentially the same basic function as the other. Air from the rotor 
enters the stator, which Is placed close behind It, usually within a distance 
of 10 to 20 percent of the blade chord. The rotor Imparts kinetic energy of 
rotation to the basically axial Incoming free-stream flow and also Increases 
Its potential energy In the form of a static pressure rise while passing 
through the Interblade passages. The stator converts the kinetic energy of 
rotation of the entering air Into potential energy by a further Increase of 
the static pressure so that the flow downstream of the stator Is, again, nearly 
axial . 

Current literature on the three-dimensional flow In turbomachines Is con- 
cerned mainly with the flow over one row of rotor blades In Isolation. In the 
single actuator disk model the perturbation velocity of the disk decreases 
exponentially with the distance from the disk. Qualitatively, the strong 
upwash field of the stator blades affects a substantial portion of the flow 
over the rotor blades and In turn Increases the effective Incidence of the 
rotor blades. Thus the rotor blades are closer to positive stall with the 
stator than without, when the effective Incidence of the rotor blades Is high 
In a positive sense. Similarly, when the effective Incidence of the stator 
blades Is high In a negative sense, the stator blades are closer to negative 
stall with the rotor than without. 

Thus It Is Important to consider the combination of the two rows, which 
form one stage of an axial turbomachine, In order to understand their Inter- 
ference effects and obtain a more accurate evaluation of the dynamic and 
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aerothermodynamlc behavior of the multistage compressor as a system through a 
synthesis of the Individual stage performances. 

Actuator disk models of the multistage compressor basically assume each 
blade row to be of zero axial thickness. In these models such an arrangement 
of many stages Is equivalent to flow through a series of thin actuator disks. 
Using the actuator disk model, Traupel (ref. 1) dealt with the case of an axi- 
ally symmetric, multistage machine with an Infinite number of Identical equi- 
distant stages. He Introduced an axially periodic stream function, with a 
period of one stage pitch, to describe the flow and obtained an expression for 
the radial velocity. 

Marble (ref. 2), Marble and Mlchelson (ref. 3), and Rail ly (ref. 4) con- 
sidered extensions of the actuator disk concept to disks of nonzero axial 
thickness. Thus Marble considered an axially symmetric flow through an actu- 
ator disk with the vortlclty shed from each blade row distributed continuously 
over the region behind the blade. He obtained a linear equation for the radial 
velocity on each side of the disk and obtained a solution for a blade row of 
finite chord by superposition. 

Rail 1 y (ref. 4) assumed the radial velocity field In the multistage com- 
pressor to be the sum of the radial velocity contribution of each stage In 
Isolation and the axial and whirl velocities to remain the same for each blade 
row. He assumed an exponential variation of the radial velocity components 
along the axis and determined the radial variation. Assuming an Initial axial 
velocity for each stage, Ra Illy calculated the radial velocity fields by super- 
position of the stage contributions and used these fields to recalculate the 
axial velocities Iteratively In order to obtain the final solution. He also 
calculated the whirl components for each case from the velocity triangles. 

Horlock (ref. 5) used the actuator disk model to study the effect of 
locating the actuator disk In the plane of the blade trailing edge and In the 
plane of the center of pressure of each blade row. Horlock and Deverson 
(ref. 6) found that theory and experiment agreed best for placement of the 
actuator disk at the mldaxlal plane of the blades. 

Kemp and Sears (refs. 7 and 8) studied the aerodynamic Interference between 
the rotor and stator-blade rows for Incompressible, nonvlscous fluids by 
regarding each blade row as an Infinite two-dimensional cascade. They obtained 
expressions for the unsteady components of lift and moment of the blades of 
each row. They also calculated the effects of stator wakes on the unsteady 
lift of rotor blades. They found lift fluctuation amplitudes of about 18 per- 
cent of the steady lift. Besides, viscous Interaction on the forces and 
moments caused unsteady forces and moments of about the same order as the 
aerodynamic Interference between the blade rows. 

Prandtl and Betz (ref. 9) outlined the lifting-line theory of the propeller 
for minimum energy losses In an Incompressible 1 nvl sc 1 d fluid. This theory 
was followed by Goldstein (ref. 10), who formulated the Incompressible poten- 
tial vortex theory of a propeller with a bound vortex line for each blade and 
a helical trailing vortex sheet shed from Its trailing edge. This theory was 
Improved and later extended to linearized compressible potential flow for the 
propeller by Busemann (ref. 11) and Davidson (ref. 12) and to the flow through 
a compressor by Rott (ref. 13). 
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A linearized three-dimensional lifting-line theory for an axial compressor 
blade row In an Infinite axial duct was proposed by McCune (refs. 14 and 15) 
and by McCune and Okorounmu (ref. 16) for both subsonic and supersonic flow. 
However, McCune's results are applicable to nonllftlng blades. Later, Namba 
(ref. 17) proposed the lifting-surface theory for a rotating thin blade row 
for subsonic and supersonic Mach numbers that uses a distribution of oscillat- 
ing pressure dipoles on the blades. His theory does not consider the effects 
of blade thickness, camber, and Incidence. The effects of trailing vortices 
shed by the blades are omitted since he used the acceleration potential. Fur- 
thermore, Namba' s theory does not consider the radial velocity or the swirl 
velocity component at the Inlet other than the circumferential velocity due to 
blade rotation. 

Wu (ref. 18) proposed a linearized, three-dimensional, compressible fluid 
flow for axial-, radial-, and mixed-flow turbomachines and outlined a numerical 
solution technique for the differential equations. 

The present report deals with the direct turbomachine problem by consider- 
ing the rotatlonally symmetric, three-dimensional, steady, Incompressible Ideal 
fluid flow through an axial compressor stage consisting of a finite number of 
blades In the rotor and stator. The rotor and stator are assumed to be located 
centrally In an Infinite, coaxial, cylindrical duct with only a small clearance 
between the blade tips and the duct walls. 

The stator experiences a periodic flow when cutting through the multiple- 
start helical vortex sheets of the rotor wake. For simplicity. It Is assumed 
that the discrete multiple-start, helical, trailing vortex sheets may be 
replaced by an equivalent continuous vortex cylinder of the same root and tip 
diameter but with uniform vortlclty over Its cross section. In this represen- 
tation the stator blades will experience a steady Incoming flow. The effect 
of nonunlform/dlscontlnuous wake vortlclty Is considered separately. 

In the present study both the rotor and the stator blades are considered 
to be straight, rigid, and untapered. The Incoming flow Into the rotor Is 
assumed to be uniform and axial with no radial or swirl component other than 
that due to the rotor rotation. Furthermore the Inflow Into the stator Is 
assumed to be primarily a uniform axial velocity with a varying swirl component 
Imparted by the rotor. The radial Inlet velocity component due to the rotor Is 
neglected at this stage. 

The undisturbed free-stream velocity components are taken to be (0, V r , W r ) 
In the (r, e, Z) directions relative to the rotor In a cylindrical coordinate 
system. The air Is assumed to enter the rotor with uniform upstream static 
pressure p^, and density p m and with an axial velocity W r that Is 
uniformly distributed over the rotor face. 

To simplify the mathematical treatment and the application of surface 
boundary conditions, the stator Is considered to be stationary. For the stator 
the radial velocity component In the rotor outflow Is neglected as negligible 
so that the Inlet velocity components for the stator are assumed to be 0, V s , 
and W s . Since the stator Is situated close to the rotor exit In a region of 
rapid change. It Is not possible to define, a priori, the Inlet velocity com- 
ponents 0, V s , and W s exactly. However, the stator Inlet conditions are 
assumed to correspond approximately to the value obtained from the velocity 
vector diagram (fig. 1) so that the Inlet velocity to the stator would have 
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the components 0, V r - W r tan <* 2 , and W r . Assuming that the rotor and stator 
are lightly loaded, this Is tantamount to the hypothesis that the perturbation 
velocity components due to the rotor are small as compared with components 
given at this stage, that the stator axial velocity W s Is uniform and equals 
W r , and that the circumferential velocity V s Is uniform and equals 
V r - W r tan a 2 . 


In the following sections a scheme for representing the lifting-surface of 
rotor and stator blades of arbitrary geometry through a distribution of flow 
singularities Is discussed, and their Induced velocity fields at an arbitrary 
point of the flow are obtained. The rotor-stator Interference factor Is Intro- 
duced next, and matching of the resultant flow field of the stage to provide 
zero net vortlclty downstream of the stage Is discussed. The boundary condi- 
tions on the blade surfaces are given In terms of blade and cascade geometry. 
These are reduced to a set of simultaneous algebraic equations to determine 
both the set of constants giving the distribution of flow singularities and 
the Interaction factor. The problem Is then discretized, and an Iterative 
scheme for the solution of the matrix of unknown constants Is outlined. The 
net pressure distribution on the blades Is expressed In terms of the Induced 
velocities, and applications are briefly discussed. 


SYMBOLS 

A expansion coefficients for rotor-blade chordwlse vortlclty 

distribution (eq. (11a)) 

A modified rotor chordwlse vortlclty distribution coefficients 

(eq. (27)) 

si matrix of coefficients A m (eq. (61)) 

AA column vector of constants A m , B m , C m , D m , and c 1 

B expansion constants for rotor-blade chordwlse source 

distribution (eq. (11c)) 

B modified rotor-blade chordwlse source distribution 

coefficients (eq. (27)) 

& matrix of coefficients B m (eq. (61)) 

C expansion constants for stator-blade chordwlse vortlclty 

1 distribution (eq. (lib)) 

C modified stator chordwlse vortlclty distribution 

coefficients (eq. (27)) 

Cp, rotor- and stator-blade half-chord lengths 

C* dimensionless rotor- and stator-blade half-chords 

# matrix of coefficients (eq. (61)) 
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E rm* E sm 
EE 

EF 


EMi 

r m 
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9r* 9s 
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1 m 




hj> , hj 

An 

k r , k s 
L r . L s 
M 
H* 

N* 


Poo 

Pr* Ps 
Por* Pos 


local lift coefficient of rotor and stator blades 

expansion coefficients for stator-blade chordwlse source 
distribution (eq. (lid)) 

modified stator-blade chordwlse source distribution 
coefficients (eq. (27)) 

matrix of coefficients (eq. (61)) 

rotor and stator functions (eq. (35)) 

matrix defined In (eq. (59)) 

matrix of Integrals (eq. (57)) 

matrix elements defined In eq. (60) (1 = 1,2, 3, 4) 

vector Integral defined In eqs. (28a) and (52) 

matrix defined In (eq. (50)) 

vector Integral defined In eqs. (28b) and (52) 

matrix defined In eq. (50) 

rotor and stator geometry functions (eq. (49)) 

vector Integral defined In eqs. (28c) and (52) 

matrix defined In eq. (50) 

rotor and stator (hub/tip) radius ratio 

vector Integral defined In eqs. (28d) and (52) 

matrix defined In eq. (50) 

rotor and stator circumferential mode numbers 

lift per unit span of rotor and stator blades 

Mach number 

matrix dimension parameter, R*N* 
number of chordwlse stations on blade 
free-stream static pressure 
perturbation pressures due to rotor and stator 
total pressures ahead of rotor and stator 
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Qr 
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R* 
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r l 

r hr» r hs 
r tr» r ts 

$rh> S sh 
SR 

T rm> T sm 
®R. ®S 

Or. Os 

-> 

u 

A 

u 

V r . V s 

^r» % 

t i 
V r » V s 


net pressure difference on rotor and stator blades 
(eq. (68)) 

total source density on rotor blade (eq. (11)) 
total source density on stator-blade (eq. (11)) 
functions defined In eq. (30) 

ratio of stator tip radius to rotor tip radius, rt$/rtr 

number of stations along blade radius 

position vector of arbitrary point (cyllndrlcal/Carteslan 
coordinates) 

dimensionless radial coordinate, r/r^ r 
hub radius of rotor and stator 
tip radius of rotor and stator 

functions defined at rotor and stator hub (eq. (36)) 
matrix defined In eq. (51) 
rotor and stator functions defined In eq. (47) 
resultant total velocity of rotor and stator 

modified resultant total velocity of rotor and stator 
(eq. (39)) 

resultant Induced velocity vector (eq. (25)) 

modified resultant Induced velocity (eqs. (26) and (27)), 
u/W a 

circumferential velocity of air for rotor and stator 
free-stream velocity vector for rotor and stator 

modified free-stream velocity vector for rotor and stator 
(eq. (39)) 


W r , W s axial velocity of air for rotor and stator 

(x,y,z) Cartesian coordinates of arbitrary point 

(y ' * z'), (y',z‘) Cartesian coordinates of rotor and stator In local coordinate 
5 system of blade 
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z r* 

Zs 

3n 


3s. 

3>s 

z Cr ’ 

Z Cs 

Z Lr ’ 

Z LS 

z Tr ' 

Z Ts 

Z Ur * 

Z US 

z r0» 

z s0 

z rl * 

z r2 

z sl. 

z s2 


a r * a s 
a 2r 
r r» r s 

C 

C 1 

e 

(k) (k ) 

\ . H 

v r , V S 

p (p. 'I'. f) 


Pr» ps 



number of blades In rotor and stator 
rotor vector function (eq. (29)) 
stator vector function (eq. (29)) 

mean-line ordinate of rotor- and stator-blade profiles 

lower surface ordinate of rotor- and stator-blade profiles 

local half-thickness of rotor- and stator-blade profiles 

upper surface ordinate of rotor- and stator-blade profiles 

axial position of midrotor and midstator plane from 
reference origin 

axial coordinate of stator-blade leading and trailing edges 
axial coordinate of stator-blade leading and trailing edges 
blade angle settings of rotor and stator blades 
exit blade angle of rotor 

total vortlclty density of rotor and stator blades 
(eq. (11)) 

rotor-stator Interaction factor 
Interaction parameter (eq. (56)), 1/(1 + c) 

azimuth angle of cylindrical coordinate system 

a-th eigenvalue for mode m r , m s 

mean azimuth angle of m r -th rotor blade, m s -th 
stator blade (eq. (14)) 

position vector of a point on rotor/stator with cylindrical 
coordinates p, 4 r, c 

radial position of source/vortex on rotor and stator blades 

normalized Bessel function of rotor blade of mode number 

k and Jl-th eigenvalue 
r 


<Ps 


azimuth angle of source/vortex on rotor and stator blades 
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mean offset angle of first rotor and stator blades 


normalized Bessel function of stator blade of mode 
number k $ and l-th eigenvalue 

azimuth angle of point on m r -th rotor and stator blades 
defined In eqs. (12) and (18) 

angular velocity of rotor 

Glauert angle of blade defined In eq. (1) 


rotor, stator 
hub, tip 

unit vectors along (x, y, z) directions 
REPRESENTATION OF LIFTING SURFACE 

The blades of the stage are considered to be thin with small camber. The 
pressure distribution on the blade may be regarded as arising from a surface 
distribution of flow singularities that give rise to the given blade profile 
and Its lift. The thickness effect of the blade Is represented by a surface 
distribution of sources and Its lift distribution by a surface distribution of 
vortices. The distribution of both the source and vortex singularities varies 
radially and chordwlse over the blade section. A study of the linearized par- 
tial differential equations for the perturbation pressures of the rotor and 
stator when the fluid Is Incompressible shows that the radial variation of the 
pressure satisfies Bessel's differential equation. Therefore, to provide for 
the radial variation of the singularities, the chordwlse variation Is modulated 
by a Bessel function. 

To specify the chordwlse distribution. It Is convenient to refer to a 
locally rotated coordinate system ( Yp - Zp) (fig. 2) and to define the 
Glauert angles w r and u s , for a point on the rotor- and stator-blade chords, 
respectively, through the equations 


<p r , <p s 



M»r* M»s 
G 

<*>p f a> s 

Subscripts: 

r,s 

h,t 

1, j. k 


y’ = -C R cos a> r ; -C R < y^ < C R ; 0 < u> r < tr 1 

( 1 ) 

y^ = C $ cos o> s ; -C $ < y^ < C $ ; 0 < g> $ < ir ) 

If the midpoint of the blade chord Is used as the reference origin for each 
blade, the unit vectors In the (R, Y', Z') coordinate system are related to 
the unit vectors (R, e, Z) In the cylindrical coordinates by the 
transformation 
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(2a) 



/ R \ / 1 0 

f Y s )=( 0 s1na s 

\ Z ' / \0 -COS a 

\ S / \ S 


° os «rY°^ 

sin « r /\ Z / 

0 V R \ 

cos a s )( 6 I 
sin a s /\z/ 


(2b) 


This system of coordinates Is rotated from the plane of rotation by the angles 
a r and , respectively. The relation between the coordinates (Y Z^) 
and (Yj, Z s ) of a point on a rotor or stator blade with the corresponding 
coordinates (Y r , Z r ) and (Y s , Z s ) In the (Y - Z) system Is given by the 
equations 




COS ¥ 

r 



sin * r 


0 

-sin a 

r 

sin 

* 

r 

sin a 

r 

cos ¥ 

r 

cos 

COS a 

sin 

¥ 

-COS a 

cos V 

sin 

r 


r 

r 

r 


cos ¥ 

s 



sin ¥ s 


0 

-sin a 

sin 

¥ 

sin a 

cos ¥ 

cos 

s 


s 

s 

s 


COS a 

sin 

V 

-COS a 

cos ¥ 

sin 

S 


s 

S 

s 




rl 

■ z r0 - C R cos V 

Z r2 " Z r0 

si 

■ z so - c s cos V 

Z s2 = Z s0 


+ C R COS a r 
+ C COS a 

v s 


(3) 


(4) 


(5) 


and (z r i, z r 2 ) and (z s -|, z S 2 ) refer to the leading and trailing edges of the 
rotor and stator. The cylindrical coordinates of an arbitrary point y^, y$ 
on the rotor- and stator-blade chord at any radius r are given by 


_ p 

r-| ; <P r + « r - arctan — cos « r sin a p ; z rQ + C R cos a> r cos a r 


S 

r l ; <p s + Uj _ ar ctan — cos sin a $ ; z $0 + C $ cos w s cos 

- 1 — Yp * y s — +1 ; h r < r < 1 ; h s R sr < r < R sr 

Note that the chordwlse coordinates y ^ and y$ are normalized 
relative to the chords 2 Cr and 2Cs, respectively. 

If zc r and zq s refer to the mean camber line of the rotor- 
and stator-blade profiles and zf r and zf s refer to the half- 
thickness of the corresponding blades, we have the relations 


( 6 ) 
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z 
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Cr " z Cr* y r ); 

z 1 = z 1 

Tr z Tr 

cs ■ z cs (y ;> : 

z 1 = z 1 

Cs z Cs 


(y;> 

(y;) 


( 7 ) 


The equation of the upper and lower surface blade profiles can be obtained from 


1 = z 1 + Z 1 • 

Ur z Cr z Tr’ 

z‘ = z 1 - z* 
z Lr z Cr z Tr 

1 = z 1 + z 1 • 

Us z Cs z Ts* 

z Ls z Cs z Ts 


Following Schllchtlng (ref. 19), the chordwlse distribution of the flow sin- 
gularities Is assumed to have the form of a Glauert-BI rnbaum series (refs. 20 
and 21). Furthermore because of the three-dimensional nature of the flow, the 
Glauert expansion will be modulated by suitable functions §* and of the 
radius. These functions are obtained as sums of the eigenfunctions 
< k r > <k s > <k) (k ) 

^ • In turn, and are normalized linear combinations 

of Bessel and Neumann functions of order k and k for the 8.-th eigenvalue 
(k ) (k ) 

and of the rotor and stator 


(k r ) (k ) 

The functions and 't are normalized linear combinations of 

Bessel and Neumann functions of order k r and k s for the 8-th eigenvalue 
for the rotor and stator and can be written 


i(k). 


Y< k ,) (>.* k) h r ) - Y (k * 1> (\J k) h r ) 
J < k - 1 >(,( k > V - J < k * 1 )(x< k >., r ) 


x * Y< k » (ri x< k >) 


't 


m 


V <k -' ) (X< k >R sr h s ) - Y< k> 1 )(y( k >R„h.) 


sr s' 


(kj (k ) 

The pressure eigenfunctions and ^ 

the root and tip expressed by the equations 


(9) 


L J (k ~ 1) (x' k) R sr h s ) - J (k * ,) (\< k) R sr h s ) 

x 0< k >(X< k ) ri ) ♦ Y (k) (X< k) r 1 ) J 

satisfy the end conditions at 
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( 10 ) 


-<L $ (k) (r ) - 
dr ] *i 'V 


A ,(k) 

dr * 


(r,) 


• 0 


at r, = h„ and r, = 1 
1 r 1 


at ri = h s R sr and r, = R $r 


which represent zero radial velocity corresponding to zero radial pressure 
gradient at the hub and tip of the rotor and stator. The eigenvalues 
(k) 

for the rotor and stator are given by the roots of the corresponding 

transcendental equations obtained by satisfying the end conditions (eq. (10)) 
at the blade tip. 


Although the Inlet flow conditions for the stator are really periodic as a 
result of the perturbations Introduced by the rotor blade passing each point 
with a frequency u s = Z r Q/2ir, In this paper, It Is simply assumed that the 
stator-blade Inlet conditions are steady by putting u s = 0. 

The total surface density of source and vortlclty for the complete range 
of eigenvalues (8. = 0,1,2,..., ») and the mode numbers k r and k s (k r , 
k s = 0, ±1, ±2,..., ±«>) Is given by 


r r^r> ■ 2W a( A 0 cot 2 


L *X A i» s ’ n m “ r ) *•<'>> 

m=l / 


V?,> - 2w a( c 0 cot ir * X 

\ m=l 


C m s1n "^s ) Y * (p s } 


Q r (P r ) = 2W a 


8 0 ( cot r - 2 s,n ^ r ) . B m s1 " 


mu 


m=2 


Q S (P S ) = 2W a 


/ (*) V 

D ol cot s1n 

m=2 


mu. 


$ *(p r ) 


^*(pJ 


(11a) 


(11b) 


(lie) 


(lid) 


00 00 


(kj 


00 00 


^ - X IC'l'rl - X X^'s* 


k =-<*> a.=o 
r 


k s =-oo suO 


The expansion coefficients A m , B m , C m , and D m (m = 0, 1 , 2, : . . , ®) are 

determined by satisfying the boundary conditions on the rotor and stator blades 

simultaneously. The cylindrical coordinates (p r , C r ) and (p s , ii? s , 

f s ) of a point on the m r -th rotor blade and m s -th stator blade are 

related to the corresponding chordwlse coordinates (u r , p r ) and (u s , p s ) the 

relations 
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2irm \ / CpS 1 n a cos u 

»l» r = «P r +1 ■y- t l - arctan 1 r 


'2irm \ / C_ sin a cos co 

t s = <p s + \~z' + arctan 5 1 


02 ) 


C = Z - + C_ COS a COS co 

r rO R r r 


C = Z A - C_ COS a COS co 
s sO S s s 


The distributions assumed In equations (2), (3), or (9) for the flow 
singularities are such that the vortlclty density r(p, cp, C) becomes Infinite 
at the blade leading edge (co = 0) to give Infinite suction and vanishes at the 
trailing edge («o = ir) , satisfying the Kutta condition as In the case of a flat 
plate at Incidence for both the rotor and the stator. The assumed source den- 
sity Q(p, cp, C) vanishes at the trailing edge (co = «) and becomes Infinite 
at the leading edge (<o = 0), corresponding to the case of symmetric Joukowskl 
profile at zero Incidence. 


Induced Velocity Field of Bound Vortlclty 

Assuming that the bound vortex filaments have their axes along radial lines 
and considering a rotor-blade surface area element (dp r df r sec a r ), the vor- 
tlclty contained In the element Is r r (p r , <p r , f r )sec a r dp r df r . The Induced 

velocity du, at the point (r, e, z) In the^flow field, due to the bound 
vortex filament with a unit vector ft at 2 (p . cp , O can be obtained 
from Blot-Savart's law as r r r r r 


,4 -4 . 

(r - pj 


du = r (p , <p , c ) sec a H x - dp df 

r r'^r’ v r' s r r r . ,•> -4 .3 p r v r 


4ir|r - r p |' 


(13) 


which can be written as 


du r = { -?(z - C r ) sin » r + J (z - C r ) cos c|» r - i<r sin (e - ^ r ) | 

r r (p r , Hy, f r ) sec a r 


4vRp( r.| , t r ) 


dp r d <r 


(14) 


where t, J, and i< are the unit vectors along the (x, y, z) directions and 
the angle ip r Is defined by 


= <P r + <P r + v p v r = 2«m r /Z r ; m r = 0, 1, 2,..., (Z p - 1) (15) 

The Induced velocity du r due to all of the Z r blades and all of the circum- 
ferential orders k r at the blade element Is given by 
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4if cos a r du p = £ [4<z - c r ) sin <!» r + J(z - f r ) cos Tj/ r - l<r sin (e - * r )J 


r r (p p * t r » c r ) 
R?(r r P r ) 


dp p df r 


Introducing the vortlclty distribution from equation (2a) Into equation 
(15) and Integrating with respect to p r and ? r gives the Induced 

velocity u r due to the vortlclty on the whole blade surface. 


2ir cos a p u rr 
W 


1 z „Z -1 


= f f EE [-?(z, - ?_) sin c + }(z - c ) 

H t m i, _ — 1 > r » 


h r Z r1 m r =0 k r = -°° 


x cos t|» r - icr sin (e - 'i' r )][A 0 cot A m sin mu r 

( 17 ) 


m=l 


x **(p r ) R; 3 (r r P r )df r dp r 


Similarly, by defining <j»s and »s by 


tlj = <p + <p + V 

S 


2irm 


m s = 0, 1, 2,...m (Z $ - 1) 


(18) 


V _ = 


s " Z. 


We can write the Induced velocity due to the bound vortlclty on the stator 
blades as 

2ir cos a $ u $r rsr /s2 “ Z - _1 


W. 


-I K f Jin t4<Z l - C S> Sl " *5 * J(Z 1 - f S> 


R sr h s z sl k s=-” V° 


x cos »|» s - kr sin (e - t s )][C 0 cot 


x **<p s )R; 3 <V * s )df s dp s 


?*E 


C_ sin mu 


m 


m=l 


(19) 


Induced Velocity Field of Source Distribution 

For a rotor-blade surface-area element (dp r df r sec o r ) , the source 
of strength Q r (p r , <|i r , f r )sec a r dp r d< r at (p r , <j> r , f r ) Induces 
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I I I 

a velocity du r (u r , v r , w r ) at the point (r lt e, z) In the flow. The veloc- 
ity Is obtained from Blot-Savart's law as 


r i - p . 


du = Q (p * , C ) sec a — 3 * - dp df 

P 4irRj <r r p r ) r r 


( 20 ) 


whose Cartesian components can be obtained from the expression 


4ir cos a du„ 
r r 


[(r-jCos 0 - p r cos + (r sin e - p r sin <j» r )J + (z 

Q r (p r * 4» r * C r ) 

— r ' dp r dC r 


C r )£] 


R?(r, p r ) 


( 21 ) 


The Induced velocity du r due to the sources on all of the Z r blades located 
at the same relative location for the whole set of circumferential orders k r 
Is given by the summation 


4ir cos a du 
r r 

V 1 


m =0 
r 


f(r 1 cos e - p r cos »j/ r ) + J(r sin e - p p sin »|» r ) + ic(z ] - 


<r> 


»r<V ?r> 3 


x Q r (p r , »i» r , C r )dp r dC r 


( 22 ) 


Again, by Introducing the source distribution from equation (lie) Into 

equation (22), the Induced velocity u r due to the source distribution over 
the whole rotor-blade surface can be obtained by Integrating with respect to 
p r and Cr an d Is given by 


2ir cos a\ 
) 



^ [t(r-, cos e 
m r =0 


P r cos t r ) + J(r sin 0 


p t sin ^ r ) + i<(z 1 - f r )J 

B 0\ cot 2~ - 2 s1n u r) + Y] B m s1n % 


X m =2 

O V V 


* P r R r (r r P r )dp r df r 

(23) 


Similarly, the Induced velocity due to the source distribution on the Z r 
stator blades can be written for the whole set of circumferential orders k s 
as 
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2ir COS a 


R .sr Ul V 1 


w_ 


3 


/ J I /, d 

7 [f(r 1 cos e - p s cos 4 > s ) + ftr sin e 


R sr h s Z s1 V° 


- p s sin * s ) + ic(z 1 - f $ )J 


x '?P s R s 3 (r- ) , P s )dp $ df s 


D q cot j- - 2 sin w$ + D m sin mu> s 


m=2 


(24) 


Induced Velocity of Combined Source and Vortex System 
of Rotor and Stator 

For low subsonic axial flow (H « 1) the resultant Induced velocity u (r) 
at an arbitrary point r (r, 0, z) of the flow field due to the combined 
system of sources and vortices on both the rotor and stator Is obtained from 
the equation 


u = 3 r r + “sr + *Vq + “sQ 


(25) 


By combining equations (16), (18), (23), and (24), u can be expressed In terms 
of the expansion coefficients of the Glauert series of equations (2) as 


u = V (A? + b3 + C 3 + D J ) 

/ v v m m mm ''mm mm' 


m=0 

where u, A m , B m , C m , and D m are redefined by the coefficients 


A„ = A /cos a 
mm r 

C m = C /cos a 
mm s 


B m = B /cos a 
mm r 

B = D /cos cl 
mm s 


u = u/W. 


Q = Q/2W. 


and ? m , ^m* ^m* an d ?m are the set of vector functions defined for 
the rotor and stator by 


(26) 
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= = === |«(P r )<lo r 'l { r 

sin mu 


0) 

r 

c °t 2~ 


4 .■* •> 


“r 

cot 2~ 


r i * p r ^ l = = === ) ® ( Pp ) dp r dC r 
sin mu 


(28a) 


(28b) 


Similarly, for the stator the vector functions are 



Mr 


. . .cot — , 

W ?s> ===== n(p s )dp s dc s 

sin mu 

Si 


Up 

cot ^ 


*2s<V V = = === *. ( P s )d 's d ‘s 

sin mu $ 


-i(z - C r )s1n t r + ?(z - f r )cos - icr s1n(e - ^ r ) 


f lr 


e 2r 


«A- ?r> 


t(r cos e - p r cos t r ) + J(r sin e - p r sin ^ r ) + i<(z - c p ) 


R ?<V p r> 


_> -t(z - f s )s1n i|» s + l(z - C s )cos - £r s1n(e - »|» s ) 

^ ls R 3 fr o 1 

R s (r l’ V 

t(r cos e - p s cos * ) + l(r sin 6 - p sin y ) + i<(z - f ) 


*2s 

with 


R s<V *s> 


(28c) 


( 28d) 


(29) 
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R r (r r p r ) = [(r 1 

R s<V p s ) - [<r, 

Besides considering the Induced velocity of the bound vortices In the rotor 
and stator, It Is also necessary to Include the Induced velocity of the trail- 
ing vortices shed by the rotor and stator. These vortices move downstream 
along helical paths. The trailing vortices shed by the stator blades are 
opposite In their sense of rotation to those shed by the rotor blades. Fur- 
thermore, when the rotor and stator have nearly equal reactions, the trailing 
vortices of both the rotor and stator may be assumed to be nearly equal In 
strength at all radii. Thus the net vortlclty downstream of the stator Is 
assumed to be nearly zero. Furthermore, because of the close spacing between 
the rotor and the stator, the effect of the vortices shed by the rotor In the 
rotor-stator gap will also be neglected. Consequently the 

4 

Induced velocity u of equation (26) Is the complete Induced velocity of the 
stage. 


- P r ) 2 + ( - C r ) 2 - 4r Pp cos 2 ^(6 - i|r r )] 

- p s ) 2 + (z-, - C s ) 2 - 4rp s cos 2 |(e - v $ )] 


1/2 


1/2 


( 30 ) 


MATCHING OF ROTOR AND STATOR FLOW FIELDS 


The perturbation velocities u (r, t) and the resulting perturbation 

pressures given In the preceding section at a point r In the flow field are 
based on the hypothesis that the Individual contributions of the rotor and 
stator are additive and are not vitiated by the Interaction effects between 
them. The singularity distributions on the blade surfaces assumed In 
equations (2) Imply no rotor-stator Interference and satisfy. Individually, 
the Kutta condition at the trailing edge of both the rotor and stator. How- 
ever, the flow tangency condition Is disturbed when the rotor and stator are 
juxtaposed closely to form a compressor stage. 

It Is necessary to satisfy the flow tangency condition for the combination 
simultaneously. This can be done by Introducing an Interference velocity 
-» 

u rs and the corresponding Interference pressure p rs so that the 

resultant Induced velocity vector v and the resultant perturbation pressure, 
p are written as 

■» -» 

AAA A . . _ . 

v = u + u rs ; p = P + P r$ (31) 


Under the hypothesis of small perturbations the Interference velocity u rs 

-4 

Is assumed to be a constant fraction c of the Induced velocity u of the 
rotor and stator 


V = (1 + e)u (32) 
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From the lifting-line model It Is known that the strength of the bound vor- 
tex Is a maximum at the blade root for both the rotor and the stator. A vortex 
filament of this maximum strength extends downstream to Infinity from the root 
of each blade. If r r h and r S h be the strength of the blade root 
vortlclty per blade of the rotor and stator, respectively, the condition of 
zero net vortlclty behind the stator blade becomes 


Vrh 


* V.b 


(33) 


This equation provides the necessary condition for determining the set of 
expansion coefficients In the Glauert series and the Interaction factor e. 
The vortex strengths r r ^ and r S h are obtained by a chordwlse Integration 
of the surface density of vortlclty In equation (11) so that we have 


r . = £ A EJhJ; r . = X C m E cm (h c ) 
rh ±—‘ n m rm' r' sh j~ n m sm s' 

m=0 m=0 


where E rm (h r ) and E sm (h s ) are defined by 


E rO< h r> " S 


rh 


E sO< h s> " S sh 


ir 

i 

i 


cot ^-dc; E (hj = S 


<0„ 


rm' r 


rh 


cot— df c ; E sm (h s ) = S 


sh 


ir 

i 

i 


sin mu, dr 
r r 


sin ma> s dc $ 


(34) 


(35) 


with S r h and S s h obtained from equation (9) by setting r-| equal to 
h r and R sr h s for the rotor and stator so that 

S . (h ) = $.(h ); S . (R h ) = (R h ) 
rh v r *' r sh v sr s' * v sr s' 

Equation (35) can be Integrated and rewritten as 

E rO* h r* = ~ 2irC R S rh* V C0S a r ; E sO* h s* = _2irC s S sh^ R sr h s^ cos a s 

E rl <h r> ■ -wC R S rh^ h r )cos V E sl (h s> = -* C s S sh (R sr h s )cos “s 

E rm<V ■ 0; E S1 »<V * 0 m 3 - 


(36) 


(37) 


BOUNDARY CONDITIONS 


For compressors and fans with small tip clearances, the combined perturba- 
*) 

tlon velocity field v of equation (32) together with the gross velocities of 
the free stream must satisfy the condition of no radial flow at both the hub 
and tip of the rotor and stator. Furthermore, both the rotor- and stator-blade 
surfaces must be stream surfaces. The latter condition Is convenient to apply 
while dealing with the resultant velocity field of the flow singularities. 
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The gross free-stream velocities \^r (0, Vr, W r ) and (0» Vs, W$) 
for the rotor and stator assumed the Introduction together with the 

resultant perturbation velocity v of equation (32) give the resultant total 

velocity Or, 0 s 


Or = 9r + v; 0 S = V s + v (38) 


where 0 R , 6 S , 0 R , and V s 
that 


are velocities renormalized In such a way 


\ K ■ 9 s /W a 


^ D /w, ; V c = t/w 


'R " R a 


s a 


(39) 


whose Cartesian components for any azimuth angle 4 » r , tj» s on the rotor 
and stator are given by 


1. - 57 s1n V - t C0S V 1 


v s ■ 


1 


1 


— + tan a 0 1 sin * e ; I r - + tan a ? I cos 1 


R 


U R - 


+ 

P 1 , 

r; s1n *r + V X ; 


• -1 

R 1 cos + v y ; 1 + v z 


U S " 


1 


1 


(40) 


— + tan a 0 ) sin + v v ; ( r- + tan a, cos i|» c + v v ; 1 + v ; 


R. 


where R + = (W a /ftrt r ) Is a characteristic dimensionless radius. R + may also 
be regarded as the advance ratio of the blade based on the tip radius of the 
rotor. Assuming that the Incoming flow Is uniform In front of the rotor, the 
flow characteristics at any point of the blade surface are only peculiar to Its 
radial location on an arbitrary blade. The same condition applies,, likewise, 
to the stator blade. Therefore the angles tr an< ^ M's given by equations 
(14) and (17) can be replaced by <p r and <p s , respectively, and the midchord 
line of the blade may be considered to be parallel to the X-axis. 


-> 4 

The resultant velocities u R and u s can be resolved along and 
perpendicular to the blade chord In terms of the local coordinates (Y 1 - Z') 
mentioned In the section Representation of Lifting Surface. To be consistent 
with the postulates of thln-aerofoll theory, the boundary condition Is applied 
at the blade chord. The resultant velocities are given by 
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Ry 


. - (^- * v x sin <p r + v y cos <p r j sin o p + (1 + v z )cos <x f 
Rz . - (^- + v x s1n <p r + v y cos <p r jcos <* r - (1 + v z )s1n <* r 


1 


u $y ' = - I + tan o 2 + v x s1n <p $ + v y cos <p c ) sin <* c + (1 + vjcos a c 


( 41 ) 


°Sz' = ( R 1 + tan 


a 2 + V 1n v s + ^y cos * s cos a s + (1 + v z )s1n a $ 


wherein we have neglected terms of second order and higher In v x , Vy, v z . 

The resultant perturbation velocity components vr and v s at the rotor 
and stator In the local coordinate system are related to those In the Cartesian 
system by the matrix transformation 


1 

s*. 

<> 


COS <p 

r 


sin 

^r 


o 


1 

<> 

X 

1 

v Ry ' 


sin <p r 

sin a r 

cos 

*r 

sin a p 

COS a p 



°Rz' 

= 

sin <p r 

COS a r 

cos 


COS a r 

-sin a r 


* 

v y 

v Sr 


COS q> s 


sin 

*s 


0 



A. 

V Cwl 


-sin <0 

sin a 

-cos 


sin a 

COS a 


A 

V 

Sy' 


s 

s 


s 

s 

S 


z 

a 

V Sz‘ 


sin <p 

s 

COS a 

S 

cos 

^s 

COS a 

S 

Sin a 

S 












— -/ 

ig the 

perturbation veloc 

Itles 

VR 

and 

v s In the blade 


(42) 


coordinates at the rotor and stator given by the column matrix on the left of 
equation (42), It Is to be noted that the Cartesian velocity components (v x , 
Vy, v z ) are evaluated at the corresponding points on the rotor and stator, 
respectively. The klnematlcal flow conditions of flow tangency on the blade 
surface at any radius r-j can now be expressed as 


'dz n 
c 

dy 1 


dz ’ 
c 

,dy', 


u 


Rz' 


Ry* 


Sz' 


Sy' 


z'=0 

r 


z s =0 


.dy 7 


'dz^ 

dy 7 , 


1 _ 
2 * 
C U 


Ry' 


i \ 


'Ry' 


z r’° 


z s’° 


(43) 


To be consistent with the postulates of thln-aerofoll theory. In equations (43) 
the boundary conditions are satisfied at the local chord of the blade by set- 
ting z|. = 0 or Zj = 0 as appropriate to the blade In question. 
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From these considerations and equation (12) at any given radius a relation 
between the azimuth angles (<p r , z r , r-|) and (<p s , z s , r-|) Is obtained as 


(Z r0 - z r )tan a r <* s - z s 0 )tan “s 

tan <p p = ; tan = — 


(44) 


1 '1 
These can be substituted for z r and z s In the expressions of equations (43). 

The slopes of the mean camber line and the thickness distribution for the 
rotor and stator blades as given by the derivatives on the left of equa- 
tions (43) are known for given blade profiles. These derivatives are denoted 
as 


dz ; 


dz 1 ' 
c 


OZj 1 


rc ” \dy 7 r ’ T rT ” 2 \dy ' J r 


sc “ \dy' i ’ t sT ~ 2 \dy 


<dz| 


(45) 


Combining equations (42) and (43) with equation (3) gives, after some 
simplification. 


sin a + t COS a ' 
r rc r 


if ♦ v x s,n *r * v y cos *r - <! + v zHco 5 « r - \r *<" « r j " ° 

\ r rc r / 

r 1 /sin a - t cos a \ 

r * tan “2 * 5 X s1n *s * 5 , cos *s + (1 * M cos . t” sin J s )- 0 


v x $ ' n *r * v y cos 


+ tan <*2 + v x s1n <p $ + v y cos <p s )tan <p s + (1 + v z ) = - 


<p )tan <p + (1 + v ) = Y, B m T (r,) 
r/ r z m rm l 

m=0 


Z Vm< r l> 

m=u 


(46) 


where 


t(r l> ... “r. 


*< r l> ... “s 


T rO (l V " t.t COt 2 ’ T S 0< r l> " t T cot 2 


T rT 

*(r,) 


sT 

nr,) 


T( r, ) = sin mu : T fm (r, ) = sin mu 

rm 1 -r - r sm r t_ t s 


T rl " " 2 


l rT 
i(r^) 


sT 

¥( r-j ) 


sin u • T , = -2 — sin u, 

Vt r sl t sT s 


( 47 ) 
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DETERMINATION OF CONSTANTS 


The nature of the distribution of the flow singularities on the blade sur- 


face depends on the coefficients A m , B m , C m , and D m (m = 0, 1, 2, 


») In 


equations (2). To determine these coefficients, the surface boundary condi- 
tions In equations (46) have been Introduced. Substituting for the resultant 
Induced velocities v x , 0 y , and 0 Z from equations (32) and (26) and combining 
the different quantities give for each point on the blade the following set of 
simultaneous algebraic equations: 


00 

r /R g 

23 VTm + K^\m + V*?m + V^m + 1 + c = 0 


m=0 


L * *U m * 

m=0 


r 1 /R + + tan <* 2 - g $ 

1 + c 


= 0 


(48) 


1 + r./R.tan cr 
l + r 


Z B T 
m rm 

. 1 + c 
iu= u 

oo A 

0 T 

> m sm 
= “ f-i 1 + e 


m=0 


A £ + B + D J. 

f-i m 4m m 4m nf 4m m^m 


m=0 


'3m 

+ - 

1 + c 

1 

+ 

r^/R + + tan 

+ 


1 + c 

the 

relations 

sin 

a 

+ t COS a 

r 

rc r 

cos 

a 

- t sin a 


r 

rc r 

sin 

a 

+ t COS a 


S 

sc s 


m=0 


I 


Q (co ) = 

S' COS ot — x 


rr . sin a 
rc r 


(49) 


and the set of functions 2F, and J by the matrix relations 


.F = SR * F ; S? = SR * G 
m m m m 


JtL = Sp * H ; 
m m 


r m - sr * j 

m m 


(50) 


The matrix SR and the column vectors F m , G m , H m , J m , and 

Sm are given by 


SR = 


sin 


<p. 


sin «p s 

sin <p tan a 
r r 

sin q> tan a 
s s 


cos <p r -g r 

cos <p s +9 S 

cos (p tan a 1 

r r 

cos <p tan a -1 

s s 


(51) 
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.^4m_ 
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.<Am. 


(52) 


(53) 


for m = 0, 1, 2,..., oo. The coefficients A m , B m , C m , and D m can be deter- 
mined from the set of equations (48). However, since the Interaction parameter 
c Is unknown, an additional equation must be provided. This Is done by 
adjoining equation (33) to equation (48), which provides the matching of the 
rotor and the stator. By combining equations (33) and (34), the additional 
equation becomes 


£ I" AZE + C Z E 1 = 0 (54 

„ n I m r rm ms sml 
m=0 L J 

with E rm (h r ) and E sm (h s ) defined by equations (35). Thus equations (48) 
and (54) are the set of slmultaneous^equatlorjs to be solved In order to deter- 
mine the Infinite set of constants A m , B m , C m , and D m and the Interference 
parameter e. 


Discretization of Problem 

It Is segn from equations (48) and (54) that the Infinite set of coeffi- 
cients A m , B m , C m , and Is, Indeed, required to describe the flow over 

the rotor- and stator-blade surfaces for each value of ri and results In 
Infinite matrices. The dimensions of the problem can be reduced by selecting a 
finite number of terms m = (M* - 1) In the Glauert series of equations (2) for 
both the rotor and the stator so that 4M* + 1 unknown constants must be 
determined. As for the blade surface, R* stations over the blade length and 
N* points along the blade chord are considered at which the camber and th 1 ck- 
ness profiles and their slopes are specified for both the rotor and stator. 

Thus the number of terms In each of equations (2) equals H* = R*N* and the 
matrix of the coefficients ft, B, C, and t) m Is of order 4M*. 

The chordwlse location of the points can be obtained by using the 3/4-chord 
theorem for each chord segment. Thus, since the midpoint of the chord has been 
chosen for each chord reference for N* = 4, these points will be located at 
(-5/16, -1/16, 3/16, 7/16) of the dimensionless rotor and stator chords Cr and 
C$, respectively. For R* = 3, the radial location of the points would be 
ri = [h r , (1 + h r )/2, 1] on the rotor and ri = [Rs r h s * R srO + h s )/2, R sr ] on 
the stator, corresponding to the hub, mean, and tip, respectively, of each 
blade. 
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The set of equations (48) and (54) can be written as a single matrix 
equation 


EF * AA = (EE * AA) c] 


( 55 ) 


where EF Is a (4M* + l)-order square matrix of the Integrals,^ & and J 
and the blade geometry parameters; AA Is the (4M* + l)-order column vector of 
the constants 8 m , § m , C m , f) m , and e -\ , where 


C 1 “ 1 + c 


(56) 


and the matrices EF, AA, and EE are defined by 


EF = 


EE = 


r 

9 


/ | 

E 

0 

Z E 


0 ' 

m rm 


m 

sm 
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J* 

9 

<<? 

9 


0 

0 

0 
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0 

0 

0 

^rm 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

T 0 

sm 

0 

0 

0 

0 

0 


! EM 


(57) 


(58) 


(59) 


The submafrlces F, Jf, and J are each of order (4M* x M*),; EH Is a 
(4M* x l)-order column matrix with each block of elements EM-| , EM 2 , EM 3 , and 
EM 4 given by 


EM 1 - (£ - *r)s EM 2 - (S7 


+ tan a 2r + g $ 


(60) 


EM 3 -( 

r l 

1 + R 
v R + 

tan « r ^; 

» 4 - 1 * 1 

* 

* 

and 


_ a 

A 

A 

A 

stf = 

[A q . . 

• V-i ]; 

[Bq. 

• • V-l 


_ A 

A 

A 

A 

<<? = 

C c 0 - •- 

' C H*-1 ]; 

[Dp- 

* • °M*-1 

are (1 x H*)-order 

row vectors; and 

Z r E rm and 

Zs^sm ' 

vectors. The matrix EE 

has only two nonnull 

submatrlci 


2 r 


tanj< 


(61) 


~ T sn/ T sT» each of order (M* x M*). 

It Is observed from equations (48) that the unknown Interaction parameter 
e 1 multiplies the unknown coefficient matrix AA of which c-j Is also an 
element. In this sense, the matrix equation (55) Is nonlinear. The equation 
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can be solved by an Iteration process that assumes e < 1 and expands in 
powers of e so that 


= ( 1 + c) ^ = 1 - g + + ... (62) 

Equation (55) can also be written 

EF * AA = (EE * AA) ( 1 - c + c 2 - c 3 + ...) (63) 

Bv setting c = 0, the zeroth-order solution AA(0) Is obtained from the 
eigenvectors AA of the matrix equation 

EF * AA = EE * AA (64) 

from which to zeroth order is obtained. This can be used to obtain 

the first-order ej^, and a first-order vector AAO) can be obtained 
from the equation 

EF * AA = (1 - c (1) )(EE * AA) (65) 

This matrix Iteration process can be continued until the changes In the suc- 
ceeding values of e-| or e and the eigenvectors are within acceptable 
limits. The eigenvectors provide the constants used In the distribution of 
flow singularities over the blade surface. 


NET PRESSURE DISTRIBUTION ON LIFTING SURFACE 

From the results obtained In the preceding section for u and v, Includ- 
ing the effects of rotor-stator Interaction, It Is possible to obtain the local 
static pressure on the blade. Thus, If pg r and pg s be the total pres- 
sures ahead of the rotor and stator with the corresponding air density p, the 
Bernoulli equation gives 


“Or - "r = »0s ‘ ' «,) < 66 > 

where V r and V s are the corresponding local velocities. Using equations ' 
(40) for the resultant velocities Ur and Us at the rotor and stator and 
neglecting the quadratic terms v£, v£, and v£ give 



The net pressure distribution on the blades, defined as the difference 
between the upper and lower surfaces of the blades, can be obtained from 
equations (67) to the first order In the Induced velocities v x , v y , and v z 
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T 1 ’ (i<7 s1n “r * cos “r)(^yL ‘ # Ju)- (£ cos “r ‘ s '" “r)(*k " 5 k) 


‘P. 


- tan a 0 ) sin a + COS a, 
R + 2 s s 


! - tan a 0 

1 "+ d > 


"i L - V 


COS a, - sin a 
S S 


( 68 ) 


< 5 k - 5 k> 


where vy and vl are, respectively, the total Induced-velocity vectors on 
the upper and lower surfaces of the blades evaluated at the chordllne of the 
rotor and stator blades as required. 

The local lift coefficients C lr and C ls of the rotor and stator 
blades are defined by 


ir 2 2 

q(l + r 1 VR;)2C r 


lS Q(1 + rJ/R^)2C s 


where L r and L s are the local lift per unit span of the rotor and 
stator blades and can be expressed In terms of Ap r and Ap s as 
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(69) 


(70) 


Since the flow field of the stage Is complex. It would be convenient to 
define the upwash as the axial component of the Induced velocity. Because of 
the nature of the chordwlse distribution of circulation given In equa- 
tions (9), the magnitude of the upwash velocity on the rotor depends on the 
chordwlse position of the rotor point considered. Let us consider the upwash 
velocity at the midpoint of the rotor-blade chord at the median plane. From 
equation (26) the upwash velocity Is obtained as 


u 7 = V ( AF + BG + CH + DJ) I (71) 

z / v m mz m mz m mz m mz | v 7 

m=0 1 Z=Z R 

which Is a function of the radial position r-| along the blade. Since the 
blade loading Increases toward the blade tip, the blade tips will probably be 
closer to stall with the stator than without. 


DISCUSSION 

This report Is primarily of a theoretical nature, outlining the method- 
ology for Including the stator of a turbomachine to make a combined study of 
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the rotor and stator as a subsystem. The lifting-surface theory outlined here 
provides a proper framework for the analysis. The theory can be applied to 
several Interesting cases. Thus the case of an Isolated rotor can be discussed 
by putting Z s = 0. The solution for a single actuator disk can be obtained by 
letting Z r •* ® and Z s = 0 while allowing the lift force per rotor blade to 
tend toward zero. The flow field of a pair of Infinite, two-dimensional cas- 
cades In parallel Is obtained for h r -» 1 ; h s -» 1 . 

It Is seen from the method used to represent the lifting-surface of the 
rotor and stator that each additional row of blades Introduces two more sets of 
coefficients In the corresponding Glauert series expansion. The overall aero- 
dynamic Interaction effect of additional rows on the first row can still be 
represented by e. Keeping the same number of R* stations over the blade 
length and N* points along the blade chord for specifying the blade surface 
geometry results In the size of the matrix Involved In determining the Glauert 
series coefficients being 2 x number of rows x R*N* + 1. 

The Bessel functions employed above In the distribution of the flow singu- 
larities extend to very high orders, for which asymptotic representations are 
Important for numerical evaluation. This aspect will be discussed along with 
the results for a stage of given geometry and flow condition and compared with 
measurements In a separate report. 


CONCLUDING REMARKS 

The application of the lifting-surface theory for a complete stage of a 
turbomachine of arbitrary camber, thickness, and other cascade geometry param- 
eters has been demonstrated for arbitrary flow conditions with subsonic axial 
flow. The separation of the rotor-stator Interference effect has also been 
shown. Expressions have been given for the spanwlse loading on the Individual 
blades for uniform steady Inlet flow. 
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